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We propose a phenomenological model, comprising a microscopic SO (5) model plus the on-site 
Hubbard interaction U ("projected SO(5) model") to understand the interrelation between the d- 
wave-gap modulation observed by recent angle-resolved photoemission experiments in the insulating 
antiferromagnet Ca2Cu02Cl2 and the d-wave gap of high-T c superconducting materials. The on-site 
interaction U is important in order to produce a Mott gap of the correct order of magnitude, which 
would be absent in an exact SO(5) theory. The projected SO(5) -model explains the gap charac- 
teristics, namely both the symmetry and the different order of magnitude of the gap modulations 
between the antiferromagnetic and the superconducting phases. Furthermore, it is shown that the 
projected SO(5) theory can provide an explanation for a recent observation [E. Pavarini et al., Phys. 
Rev. Lett. 87, 47003 (2001)], i. e. that the maximum T c observed in a large variety of high-T c 
cuprates scales with the next-nearest-neighbor hopping matrix element t' . 
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I. INTRODUCTION 

The SO (5) theory has been proposed as a unified 
description of antifcrromagnetism (AF) and ekwave su- 
perconductivity (SC) in the high-T c materialsa. While 
originally formulated as a phenomenological effective 
field theory for the global phase diagram of the high-T c 
cuprates, it was soon realized that the SO (5) symmetry, 
in contrast to the more common spin SU(2) symmetry, is 
not realized exactly in strongly-correlated models, such 
as the t — J and Hubbard model. However, it can be 
detected, under some|-cijcumstances, in the low-energy 
sector of such modelsolJ. This seems to indicate that 
this symmetry, which is broken at the microscopic level, 
may be restored, at least approximately, at long distances 
and low energies. 

Recently, it was pointed out that an exact formula- 
tion of the SO (5) theory cannot properly account for the 
large Mott-insulating gap which is present at half fill- 
ing in the high-T c cuprates: an SO(5) transformation 
"rotates" spin into charge and, thus, an exactly SO (5) 
-invariant system should have the same charge and spin 
gap, in contrast to the actual situation for the high-T c 
materials. This prompted for an improved formulation 
of the original idea by introducing a Gutzwiller "projec- 
tion" onto states where the coastraint of no double occu- 
pancy is implemented locally&El. In particular, a so-called 
"projected" SO(5) bosonic model has been explicitly con- 
structed:]. It is built up out of five bosons, three triplets 
which "condense" at low temperature into the antifer- 
romagnetic state and two charged bosons, i. e. Cooper 
pairs, which condense into the superconducting state. 
The "high-energy" doubly occupied charged boson states 
are then projected out from the Hilbert space. The cen- 
tral hypothesis of the projected-SO(5) (p-SO(5) ) theory 
is that the model is accurate in describing both the static 



and dynamic properties in the high-T c superconductors 
(HTSC). 

On the basis of a numerically exact Quantum-Monte- 
Carlo (QMC) calculation of the p-SO(5) bosonic model, 
it has recently indeed been shown that this model gives 
a realistic description of the phase diagram of the HTSC 
materials and-oroperly accounts for many of their physi- 
cal propertiesE3. Among those are the neutron-scattering 
resonance, which appears as a Goldstone mode in the p- 
SO(5) description, as well as an unusual chemical poten- 
tial dependence on doping found in the prototype ma- 
terial La2- x Sr x CuO$3, wiich signals a possible (mi- 
croscopic) phase separatior£3. One crucial point of this 
bosonic p-SO(5) model is that, on the basis of a newly im- 
plemented Stochastic Series Expansion QMC technique, 
it can be simulated up to unprecedented system sizes of 
about 10 4 sites, which is more than one order of mag- 
nitude larger than corresponding fermionic QMC sim- 
ulations. This fact allowed, for the first time, for an 
approximate "finite-size-study" of the p-SO(5) bosonic 
model in the three-dimensional case, and for extracting 
the "scaling properties" close to the bicritical (AF-SC) 
point. The numerical results point to a partial asymp- 
totic restoring of the SO (5) symmetry at this critical 
point, i. e. at long distances. This is very interesting, 
since the projection destroys the symmetry at the Hamil- 
tonian level (The Hamiltonian no longer commutes with 
the SO (5) generators). While the mean-field classical p- 
SO(5) Hamiltonian of Ref. ^conserves its SO(5) invariant 
form, quantum fluctuations break the symmetry. How- 
ever, as shown in RefJ13, at finite temperature quantum 
fluctuations become less and less important, and one can 
hope that symmetry-breaking effects due to the projec- 
tion become asymptotically irrelevant in the neighbor- 
hood of a finite-temperature critical point, which is what 
the numerical results seem to find. 

Therefore, one part of our motivation to study the p- 
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SO (5) model on the fermionic level stems from the suc- 
cess of the bosonic model. Thus, in the present paper, 
we start again from the assumption that the high-T c 
materials are characterized by an underlying approxi- 
mate SO (5) symmetry. This assumption is further sup- 
plemented by the Gutzwiller "projection", i. e. the 
Hubbard- [/-term provides an important additional ingre- 
dient to the SO (5) -symmetric part of the Hanhltonian, 
which guarantees the presence of a Mott gap at half- 
filling. In this sense, we choose our Hanailtonian as the 
simplest (namely scalar SO (5) invariants) Hamiltonian, 
which possesses such a symmetry and which displays the 
correct d-wave superconducting gap at finite doping. The 
physical basis is thus the assumption of a projected SO (5) 
symmetry, and the existence of a superconducting state 
at finite doping. The t — _J model was shown by our 
previous numerical studiestj to have (nearly) projected 
SO (5) symmetry in the bosonic sector. However, the 
t — J model cannot explain the |<i|-wave AF gap mod- 
ulation in the fermionic sector, it therefore misses an 
important piece of physics our current model contains 
(unless one introduces ad-hoc values for longer-ranged 
hoppings t',t", see below). Thus, the logic of our ap- 
proach is to construct an effective model, which matches 
the low-energy experiments, rather than trying to de- 
rive an effective model from first principle. This is the 
central idea behind a Landau-type approach to strongly- 
correlated systems, and can be very useful, although not 
complete. 

This brings us to the second part of our motivation. 
It is the relation to recent angle-resolved photoemission 
data, which found evidence for a d-wave-like modulation 
of the antiferromagnetic gap, suggesting an intimate in- 
terrelation between the AF insulator and the SC with its 
d-wave gap. In a recent lettem it was shown, that the 
projected SO(5) theory correctly describes the observed 
gap characteristics. Specifically, it accounts for the order 
of magnitude difference between the AF gap modulation 
of order J « 0.2 eV) and the SC gap (~ J/10) and is 
also consistent with the d-wave gap dispersion. Thus the 
projection is crucial, since in an exactly SO (5) -invariant 
model, the d-wave SC gap would be transformed into 
a pure |d|-wave AF gap with the same amplitude and 
without a constant (s-wave) part. In Ref. |^, the on-site 
Hubbard U has been treated by a Hartree-Fock decou- 
pling, unifying the gap interrelation in the much-used 
spin-density-wave (SDW) picture for the AF and BCS 
limit for the SC. This treatment will be relegated here to 
a short appendix (Sec. |A]), introducing some notations 
and also giving a simple illustration of the crucial differ- 
ence the "projection" , i. e. the [/-term induces in the AF 
gap modulation and in its SC counterpart. In the present 
paper, we concentrate on a more controlled calculation 
for larger U (U 3> t), i.e. a slave-boson calculation. 

This paper is organized as follows: 

In Sec. 0, we start with a summary of the experimen- 
tal results observed by angular-resolved photoemission 
spectroscopy (ARPES), which have motivated our cal- 



culation. In Sec. Ill, we show that the 7r-operators of 



the SO (5) symmetry "rotate" the d-wave symmetry of 
the SC order parameter into the |d|-wave symmetry of 
the AF order parameter. This argument is general, and 
independent of the details (in particular, it survives the 
projection). More specifically then, we discuss the pos- 
sibility of introducing a generalized form factor for 
the SO (5) symmetry, which emphasizes the importance 
of longer-ranged pairings in accordance with the experi- 
mental observations. The comparison of our results with 



ARPES experiments is given in Sec. Ill C . In Sec. IV 



we compare the mean-field result with a slave-boson ap- 
proach. The treatment of the antiferromagne tic an d of 
the s uperconducting phase are given in Sec. IV A and 
IV B , respectively. The usually employed t'-(t") fitting 
procedure in Hubbard- (t — J-) like models to ARPES 
data seems, at first sight, to be completely uncorrelated 
with the experimental finding, which emphasizes the uni- 
versal role of the magnetic energy scale J, namely the 
fact that t he insulating bandwidth itself is of order J. In 
Sec. Ill B, we show that the t' — t" -hoppings can indeed 
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arise effectively from the SO(5) part of the Hamiltonian, 
i. e. from the spin interaction term. Since in the SO (5) 
description these hoppings scale with the SC strength 
(V in Eq. [l0| ), they are also qualitatively in line, from 
a rather different point of view, with a recent observa- 
tion!^ i.e. that the maximum T c observed in a variety 
of HTSC scale with the next-nearest-neighbor hopping 
matrix element t' . This empirical observation a nd its 
p-SO(5) interpretation will be discussed in Sec 
Finally, Sec. [v| gives a summary of our results. 



II. EXPERIMENTAL GUIDELINES 



Our analysis is motivated Jsy recent ARPES experi- 
ments by the Stanford groupie, which indicate an inti- 
mate interrelation between the superconducting and an- 
tiferromagnetic state. These experiments show remnants 
of a superconducting gap in the half-filled high-T c par- 
ent compound Ca2Cu02Cl2 . The ARPES gap struc- 
ture in the AF phase is given by the dots with error 
bars in Fig. 1 of Ref. 0. It displays a d-wave-like, i.e. 
| cos k x — cos k y | -like dispersion in the one-electron spec- 
tral function A(k,ui) with respect to the lowest energy 
state at wave-vector k = (tt/2,it/2). The points show 
a d-wave dispersion function of the photoemission band 
along the edge of the magnetic Brillouin zone, which, 
by assuming a symmetric inverse photoemission band, 
can be understood as a modulation of the single-particle 
Mott-Hubbard gap. The inset of this figure presents 
the ARPES data and the dispersion features in a two- 
dimensional plot: on a line drawn from the center of the 
Brillouin zone to the experimental points, the distance 
between these points to the intersection of this line with 
the AF Brillouin zone gives the value of the dispersion 
at the k-point considered. The data closely follow the 
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\cosk x — cos ky | -behavior, depicted in the d-wave full 
line. The important point is that these photoemission 
data suggest that the ci-wave-like dispersion in the insu- 
lator is also the underlying reason for the pseudogap in 
the underdoped regime: this "high-energy" pseudogap 
of the order J ~ O.leV continuously evolves out of the 
insulating feature, as documented not only by the same 
energy scale but again by the same (i-wave dispersionEZI. 
Since, on the other hand, this high-energy feature is 
closely correlated to the superconductiwt-gap as a func- 
tion of both doping and momentunJiiJU'ta, we conclude 
that a microscopic theory should be able to explain the 
interrelation between the superconducting gap and the 
AF gap modulation. In Ref. ^, we have shown that a 
SO (5) model plus U provides such an interrelation in a 
natural way. In the present paper, we want to discuss 
this issue in more detail. 



III. MICROSCOPIC MODEL 

As discussed above, an important point missing in an 
exact SO (5) theory is the correct description of the Mott- 
gap at half-filling. For this reason, it is important to ac- 
count for the strong on site Coulomb repulsion U. This 
can naturally be done by including the usual Hubbard 
on-site interaction term U "by hand" . Formally, one 
can project out doubly-occupied sites by taking U — > oo, 
eventually. This should be a good approximation since 
these states are separated by an energy scale, which is or- 
ders of magnitude higher than the interesting low-energy 
scales set by the Neel (T/v) and by the superconducting 
(T c ) transition temperatures. The idea of such a "pro- 
jected" SO(5) theory, recently formulated in a rigorous 



fashion for bosonic degrees of freedom in Ref. |[ is visu- 
alized in Fig. [|. 

As already mentioned in the Introduction, the physical 
basis for the choice of our Hamiltonian is the assumption 
that the high-T c materials are following an underlying 
SO(5) symmetry. The hamiltonian in Eqs. (|T3]-|l7|) has 
been chosen as the simplest SO(5) -symmetric hamilto- 
nian that can properly reproduce the d-wave supercon- 
ducting state of the High-T c materials in a BCS mean- 
field description. Our starting point is the supercon- 
ducting state. We then perform an SO (5) rotation on the 
operator level that introduces the magnetic part of the 
interaction [second term in Eq. [H| . The resulting hamil- 
tonian is then SO(5) invariant. This hamiltonian must be 
further supplemented with the projection (Hubbard-[/) 
term and the usual chemical potential term. 

Our aim is to describe the peculiar ARPES results of 
Ref. 16 by the simple assumption of SO (5) symmetry 
of the microscopic Hamiltonian plus the projection pre- 
scription. The first task is, thus, to build up the simplest 
fermionic lattice Hamiltonian obeying these two require- 
ments. The part of the Hamiltonian leading to the su- 
perconducting state is dictated by experiments on the SC 
gap in the cuprates. These experiments suggest a d-wave 
superconducting gap resulting in a nearest neighbor d- 
wave BCS interaction. We also allow for more general 
forms of the SC gap function, including longer-range in- 
teractions, as observed in more accurate ARPES experi- 
ments on Bi 2 Sr 2 CaCu208+5 With the SC part of the 
Hamiltonian being fixed, the AF part is then also fixed 
by the SO(5) rotation. Therefore, once the form of the 
wavefunction in the SC state is given, the SO(5) assump- 
tion allows to make a prediction about the AF state. 
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FIG. 1. Suppose one starts from some singlet statelifcil as a description of the vacuum [e.g. a resonating-valence-bond-like 
(RVB) state], and looks for excitations on top of it. In an (ideal) exact SO(5) -symmetric description (a), triplet-, hole-pair- 
and electron-pair-excitations on top of this RVB state are equivalent. This condition is obviously violated in the cuprate 
materials at half filling where spin-excitations are gapless (due to the antiferromagnetic state), while hole-/electron-pair 
excitations have a large gap of the order of the Mott gap (b). A projected SO(5) theory, however, can restore the equivalence 
between spin and charge excitations in the hole sector. This is done by lowering the chemical potential to the top of the valence 
band, and projecting out high-energy electron-pair excitations so that hole-pair excitations become gapless (c). 



To see this, let us consider the operator creating a 
Cooper pair with vanishing total momentum: 



P(k) = - C k ^C_ k 



(1) 



where c k = (ck,|, Ck,J.) is a two-component destruction 
operator for an electron with momentum k, and a a are 
the usual Pauli matrices (a — x^^. An SO (5) rotation 
is generated by the tt operatoi-BoE] 



7r « = £ 



5k Ck+Q o- a cr v c_ k , 



with the form factor 

.9k = sgn(cosk x — cosky) 



(2) 



(3) 



The 7T operator rotates the Cooper pair into the magnon 
operator 



1 



N a (k) 

via the commutation rules 



2 C k+Q a c ' 



k i 



(4) 



and 



7r Q ,A^(k) = -i 5 a ,j3 0k P(k) , 



jrt,P(k) =ig k N a (k) 



(5) 



(6) 



From Eq. (^) and Eq. (^) it is clear that a d-wave 
Cooper pair <i k P(k) with the nearest-neighbor gap dis- 



persion 



d k = (cos k x — cos k y ) 



(7) 



transforms into a -wave-shape" magnon |d k | A^(k), 
giving a corresponding "|d|-wave" gap in the AF state. 
In its original formulation, an SO (5) rotation was defined 
by the tt operator, Eq. (||) , with the "nearest-neighbor" 
form factor Eq. (|3|) . From the above discussion it is clear 
that such a formulation is only appropriate for a similarly 
nearest-neighbor superconducting gap function, Eq. 
(0) . An extended gap, as foiled in recent ARPES experi- 
ments on Bi 2 Sr2CaCu208+«5 E2l would not be transformed 
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into the appropriate AF form by this SO(5) transforma- 
tion. However, Eq. (^|) is not the only possible form of 
the form factor which can be used. In fact, one can use 
a different SO (5) symmetry transformation with a more 
general form of g^. More specifically, it can be shown 
that the ir operator Eq. satisfies the closure require- 
ment of the SO (5) group as long as the form factor has 
the properties |<7k| = 1, and <?k+Q = — .9-k- Thus, given 
a more general form of the gap dispersion <ik satisfying 
c^k+Q = — rf-kj one can choose 



,g k = sgn d k 



(8) 



A gap function satisfying these requirement can support, 
e. g., third-neighbor terms: 

dy = b(cosk x — cosfcj,) + (1 — b)(cos3k x — cos3k y ) . (9) 

From the commutators, Eq. (|5|) and Eq. (|(|), it is thus 
clear that the generalized choice Eq. (||) produces a 
corresponding AF gap function of the form \dk\. 

As explained above, our procedure starts by taking the 
simplest fermionic lattice Hamiltonian reproducing the 
appropriate d-wave state in the SC state within a BCS 
(mean-field) description. This is given by 



H SC = H km + [A At + At A] 
where the kinetic-energy term reads 

Hkin — ^ ] £ P C p,(T C P-cr ■ 



(10) 



(11) 



Here, e p is the single-particle kinetic energy, and A the 
d-wave pairing operator 



A = ^d k P(k) 



(12) 



part of the Mott-insulating gap, as present in the high- 
T c materials. For this reason, our Hamiltonian must be 
supplemented with a term suppressing double occupation 
at half filling. This can be achieved by introducing the 
usual on-site Hubbard-repulsion term 



H, 



It.' 



(15) 



where, 



— r' c 



is the particle number with spin 



a on site r. and for convenience, we have subtracted a 
constant chemical-potential term. Clearly, for a complete 
suppression one should take the U — > oo limit. However, 
here we will consider a finite U = 8 t giving a more phys- 
ical partial projection resulting in a finite Mott gap at 
half filling. It is clear that, concerning the energy scales 
of the order J we are interested in, there is very little dif- 
ference between results at U/t = 8, and U/t — oo. This 
will be confirmed by our calculation later, see Fig. ||. 
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FIG. 2. The amplitude A mo J of the AF d-wave-like mod- 
ulation as a function of the Hubbard-interaction U for fixed 
V — 0.61 t and 6 = 1. The figure plots the results of the 
slave-boson calculation of Sec. 



with P(k) given by Eq. (g). 

The second term on the right-hand side of Eq. (110) 
describes a coupling between Cooper pairs. The SO (5) 
assumption requires that, (i), a corresponding coupling 
between magnons be present as well, which leads to the 
form 



V 



m 



where 



El°k|^(k) . 



(14) 



and, (ii), that £ p +q = — £— p , which suggests the simple 
nearest-neighbor form e p = —2t(cosp x + cosp y ). 

As discussed in the Introduction, restricting to a purely 
SO(5) -symmetric Hamiltonian like Eq. ([l3]) , would pro- 
duce a gap with nodes in the AF state and no constant 



The last term in the Hamiltonian is the usual chemical 
potential term controlling the doping, i. e. 



,t 



(13) In total, our Hamiltonian has the form 



Ht, 



H, 



SO(5) 



Hu + H„ 



(16) 



(17) 



In Ref. [l3], we have shown that the chemical poten- 
tial term partially compensates for the SO (5) symmetry 
breaking of the Hubbard term Hjj (see also schematic 
Fig. |l]c) , as also shown for the t — J model in Ref. 

The interacting part of H§q^ is a special case of the 

mawifestly SO (5) invariant part introduced by Rabello et 
alH 
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H S0(5),i = H S0(5) 



a=l k.k'.q 



(18) 



where 



{ c k,T> c k,i> 5( k ) c -k+q,t> ,g(k)cL k+Qa } , (19) 



is an SO (5) spinor, and the T matrices 



ia y 



p6 _ / <7* 

1 a y 



fj x /y I z o 
cWa/* 5 



(20) 



provide the SO(5) group structure. In particular, our ex 
actly SO (5) invariant part of the Hamiltonian, Eq. (IT - 
is given by Eq. (jlj) with the simple choice 

V 



Vi(k,k',q) = 



|dk| Mk'| <5q,Q ■ 



(21) 



16 N 

Our "minimal" projected SO(5) model contains, in 
principle, five parameters, although only two of them 
are adjustable. The adjustable parameters are the in- 
teraction term V, and b [cf. Eq. (||)]. V, due to the 
SO (5) condition, controls the magnitude of both the in- 
teractions between Cooper pairs and between magnons. 
As anticipated, we fix this parameter by fitting the ex- 
perimental magnitude of the gap in the SC phase. The 
second one, b , is needed when the SC gap deviates from 
the nearest-neighbor form. This parameter can be fitted 
by the detailed form of the gap in the SC phase. The 
other three parameters, zx, t, and U/t, are not free, /x 
is fixed by the hole density, which we take from the ex- 
periments to be < n >= 0.89E3 in the SC phase, while 
for t we take the commonly accepted value t — 0.5eV. 
Several experiments and theoretical fittings agree about 
a strong-coupling value of U of the order U/t w 8. On 
the other hand, we show that our results saturate in the 
large-J7 limit, i. e. they are practically independent of U 
for U/t > 8. 

We now proceed to illustrate the mean-field treatment 
of our Hamiltonian Eq. (|l7|). The main point of this 
treatment is to (a) show that the \d\- modulation survives 
the Gutzwiller projection. This important property fur- 
ther clarifies what is meant by "projected SO (5) sym- 
metry" in the ferroionic sector, generalizing the concept 
introduced in RcfM Furthermore, we show (b) that the 
SO (5) projection introduces a difference in the magnitude 
of the superconducting and antiferromagnetic gap mod- 
ulation which is consistent with the experiment. This 
quantitative analysis is more model dependent compared 
to point (a), but within the approximations we show the 
robustness of the result. 

Details of the SDW treatment of the antiferromagnetic 
phase and of the BCS treatment of the superconducting 
phase as well as some basic notations can be found in the 
appendix, Sec. 



A. Discussion of the SDW solution 



Eq. QAq ) yields the coupled set of self-consistent equa- 
tions: 



V y-^ ,A mod |d k | + Au 



2N 



E(k) 



U ^ A mod |d k | + Au 



2N 



E(k) 



(22) 
(23) 



where E(k) is the quasiparticle energy given in Eq. All, 
and where the sum runs again over the whole Brillouin 
zone. The behavior of Au and of ratio A moc i/ A$c 
reported in Fig. |^. As one can see, the ratio between the 
AF and SC gaps, which is unity at the exactly SO (5) - 
symmetric point U/t = 0, increases with U and saturates 
at U/t w 8. 
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FIG. 3. SDW mean-field results for Au and for the ratio of 
the d-wave-like modulated parts of the gap A mo( j/Asc (both 
at half filling) as a function of U. 



We now study the gap structure of the single-particle 
spectrum along the magnetic zone boundary. Since the 
free dispersion vanishes here, the gap becomes 



A AF (k) 



2(A V - 



- E%k) 
\du\ A m 



(24) 



We, thus, obtain a c?-wave-like modulation with ampli- 
tude 2 A mo d on top of the constant Hubbard gap 2 Ay, 
as anticipated above. This gap structure is represented 
schematically in Fig. (0). 
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FIG. 4. Schematic gap structure in the AF phase along 
the magnetic zone boundary. On a line drawn from the center 
of the Brillouin zone to any point on the curve, the distance 
from this point to the intersection of the line with the anti- 
ferromagnetic Brillouin zone boundary gives the gap at the 
k-point of the intersection. On top of the constant Hubbard 
gap, which, for convenience, has been reduced in scale, there 
is a d-wave like modulation. 



B. Physical interpretation of the SO(5) -generated 
(long-range) magnetic interaction part 



It is instructive to transform the magnon term Eq. 
( |l4| ) entering the magnetic part of the interaction (last 
term on the r.h.s. of Eq. (|l3|)) into real space. This 
gives 

m = ]Tm(r), (25) 

r 

where the Neel order parameter at site r reads 

m a (r) = IJ2 u '( r - r 2) ^ 4 ° a c r2 , (26) 



and has an extended internal structure given by the 
Fourier transform w(r) of \dk\. This function is given 
by (cf. Ref. [IJ) 



w(m a. x + n a y ) 



1 + (-1)™+" 



7r 2 [(to + n) 2 — 1] [(to — n) 2 — 1] 



(27) 



We want to show how, at the mean-field level, the ex- 
tended structure of the magnetic part of the SO (5) - 
symmetric Hamiltonian leads to effective longer-ranged 
hopping terms similar to the commonly used values. 
These terms can be explicitly seen in the mean-field 
SO (5) -symmetric Hamiltonian in real space 



SO(5)l 



MF 



E 



aw(r 



y^cl a c r , a + H sc 



(28) 



Notice that the effective hopping terms only connect sites 
within the same sublattice, since w(r — r') vanishes oth- 
erwise. Furthermore, the sign of the effective hopping 
matrix elements is opposite for spin up and down, as 
well as for the two sublattices, thus yielding an identi- 
cal sign for the majority spin species on each sublattice 
i.e. for the electrons whose spin is polarized along the 
direction of the staggered magnetic field. The processes 
involving those "majority" electrons are the most impor- 
tant ones since the ratio between the number of "major" 
to "minor" electrons is 



N, 



major 



(1 + (to»/(1-(to)), 



(29) 



where 



(m) = is the mean-field staggered magnetiza- 
tion which is of the order (to) i. 

The effective second- and third-nearest neighbor 
"hopping amplitudes" generated by the interaction Eq. 
<M) are 



t' = A mod -^ « -0.08t, 



t" = A 



mod 



9tt 2 



-0.03t. 



t"/t' = -1/3. 



(30) 

(31) 
(32) 



Here we have used the self-consistent solution A mo d ~ 
0.3t for large U and the results of Sec. A 2. These ef- 
fective second- and third-nearest neighbor hopping ele- 
ments have the same sign but are somewhat smaller than 
the ones commonly introduced in order to correctly re- 
produce the Fermi-surface topology of typical high-T c 
materials (see, e. g., Rcf. |23|). Nevertheless, here, these 
"effective" parameters are obtained without any fitting 
to the experiments by the simple assumption that leads 
us to the projected SO (5) Hamiltonian Eq. (0). This 
makes clear, in particular, that the gap modulation is re- 
lated to the magnetic energy scale J, as demonstrated in 
the ARPES experiments and repeatedly stressed in the 
literatureo Eq. (^8|) also contains longer-ranged hop- 
ping processes which, however, are much smaller than t' 
and t" . These long-range terms, however, produce the 
cusp-like feature of the antiferromagnetic gap modula- 
tion which has. been observed in ARPES experiments on 
Ca 2 Cu0 2 Cl 2 Ei 

Since the hoppings Eq. ( |30| ) scale with the strength 
of the SO (5) coupling (V in Eq. |l^), and since, at least 
in the weak-coupling regime V <C 8t the superconduct- 
ing T c increases with V , one can expect that T c is an 
increasing function of the next-nearest-neighbor hopping 
t' . This fact has been observed in a recent analysis of 
Pavarini et alx3, who have plotted the maximum T c as 
a function of the ratio t' /t for different compounds. As 
explained above, our theory provides a qualitative expla- 
nation for such a behavior. This is more quantitatively 



7 



shown in Fig. ||, which plots our T c versus t' results in 
the upper part. Here, different next-neares-neighbor hop- 
pings t' have been extracted by varying the SO (5) cou- 
pling V. For a given V, and thus t', T c has then been 
extracted from the weak-coupling BCS equation, which 
is valid for V < 8t. The lower part of Fig. eLgives finally 
the T c versus t' results of Pavarini et al&3 for a large 
variety of HTSC materials. 




100K 
To max 
75K 



1 1 1 - 

Hg Ba 2 Ca 2 Cu 3 O a ' □ H ■ 

Hg Ba 2 Ca Cu 2 6 ' -O . 

TI 2 Ba 2 Ca 2 Cu 3 O 10 k u w ■ 

TI 2 Ba 2 CaCu 2 B ' q, B 

Y Ba.CUaO; n _ 

HgBa 2 Cu0 4 h — ■— 
TI 2 Ba 2 Cu 6 c ■ 

D Pb 2 Sr 2 YCu 3 OB ■ 

La 2 Ca Cu 2 6 m Q ■ - B a La Cu 0. 

□■ La Ba 2 Cu 3 7 e ■"■ 

La 2 Cu Q," ■ ■ Bi 2 Sr 2 Cu 6 c 

Pb 2 Sr 2 Cu 2 6 b 

~ 5 B Ca 2 Cu0 2 CI 2 a 
I I I 

0.2 r _ t'/t 0.3 0.4 



FIG. 5. Correlation between T c 

and the next-nearest-neighbor hopping t' (in units of t) . The 
top figure reports the results of our calculation for, whereby 
different t' are obtained by varying V, and the superconduct- 
ing T c is evaluated via the weak-coupling BCS equation. The 
doping is fixed at 15% (< n >= 0.85), U = 8t, and 6=1. For 
comparison, the bottom figure shows the results of Ref. [l^ for 
diffcnt materials. 



experimental figure displays the peak dispersion versus 
the d-wave-like function | cos k x — cos k y \ . The very 
good fit to a straight line shows that the dispersion is 
d-wave-like with a characteristic "cusp" at momentum 
k=(7r/2,7r/2). 

Concerning the energy scales that are observed in the 
experimentt3, the authors find a modulation A moc i of 
the order the magnetic exchange coupling J which is 
about one order of magnitude higher than the experi- 
mentally observed amplitude of the superconducting gap 
Asc ~ J/10. This is a second, most important point that 
would be difficult to understand in terms of an "exact" 
SO (5) symmetry, which would be expected to preserve 
the amplitude of the modulation by "rotating" from the 
SC to the AF gap. However, the result of our projected 
SO (5) calculation is that the difference by one order of 
magnitude between the two amplitudes is correctly re- 
produced, and it is due to the projection as well. This 
can be seen from Fig. (|^), displaying the ratio A mo d/Asc 
as obtained by our mean-field calculation as a function of 
the Hubbard repulsion U. We see that for large U ~ 8t 
the AF gap modulation A mo( j is approximately four times 
larger than the superconducting gap AgCj when the lat- 
ter is also calculated at half filling. An additional factor 
two emerges from the doping dependence of Asc shown 
in Fig. (^|) so that the total ratio of A mo( j/Asc becomes 
of the order of 8. 




5 10 15 20 

%doping 



25 



FIG. 6. Mean-field result for the d-wave SC gap Asc ob- 
tained by solving the gap equation (Al) as a function of dop- 
ing (1 - (n)) for V = 0.8t and 6 = 1. At a doping of » 20%, 
Asc takes the value Asc ~ 0.05t which is about the correct 
order of magnitude for the high-T c materials (Asc ~ J/10). 



C. Discussion of ARPES experiments 

The interesting point is that the structure of the AF 
gap Eq. shown in Fig. [| is a consequence of the 

SO (5) -symmetry principle combined with the effect of 
the Hubbard interaction U. This gap structure has re- 
cently been measured in the undaped (half-filled) cuprate 
parent compound Ca2Cu02ClJl3. Fig. 1 of Ref. [)] shows 
the ARPES results for the AF gap modulation on this 
material, in comparison with our result, Eq. (p3). The 



At first sight, it might seem difficult to understand 
why the introduction of U changes the order of mag- 
nitude of the SC and AF gap modulations, despite the 
fact that both are controlled by the same energy scale V. 
This can be easily understand at weak coupling V <§; t 
and for U 3> t. In this case, the usual BCS gap equa- 
tion Eq. (Al) - which does not depend on U - gives 
Asc ~ exp — [const./ (npV)], i. e. it decreases exponen- 
tially with small V. On the other hand, the AF gap 
equation, Eq. (E2h, does depend on U. In fact, for large 
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U, E(k) w Au oc U, which inserted in Eq. gives 
A mo d oc V, i. e. A moc ; behaves linearly with V, totally 
different from A$c- 

Our "projected" SO (5) symmetry principle, thus, can 
provide an explanation of the interrelation between the 
SC and the AF gap both qualitatively, via the same dis- 
persion, as well as semiquantitatively, i. e. via the correct 
order of magnitude difference between the two modula- 
tions. 



IV. SLAVE-BOSON APPROACH 

The Hartree-Fock treatment described in the previous 
section was physically motivated as a natural first step 
to show the connection of the SDW with the BCS solu- 
tions via a SO (5) theory. Our motivation for introducing 
the Hubbard repulsion U was to project out states con- 
taining doubly occupied sites and therefore to account for 
the Mott gap at half filling. However, within the Hartree- 
Fock treatment discussed above, the Hubbard interaction 
is treated only perturbatively, which might be question- 
able for large values of U. For this reason, we repeat the 
calculation using an extension of the Kotliar-Ruckenstein 
slave-boson approachEjcj which, while still requiring a 
mean-field approximation, is more appropriate to deal 
with the strong repulsion U. 

Confidence in this approximation derives from various 
observations: (i) This method is known to give rather sat- 
isfying agreement with QMC results over a wide range of 
Coulomb correlations U and values for the dopingpl. (ii) 
The two main results of our study, namely, that the AF 
gap modulation is correlated to the SC gap by symmetry 
and that the Hubbard term Hjj induces the order of mag- 
nitude difference in these effects should be independent 
of specific approximations, at least qualitatively. This 
will be explicitly verified by comparing with the simple 
Hartree-Fock study described above. 

Within the Slave-boson formalism for the Hubbard 
model one introduces four bosons er , Pr]a, and 
for each lattice site r, corresponding to empty, singly 
occupied site with spin a, and to doubly-occupied sites, 
respectively. Creation and annihilation of a fermion c ri<T 
is mapped onto the creation and annihilation of a pseudo 
fermion f Tt(T with an appropriate boson, according to the 
mapping: 

c t > ft z i 

c r.cr * z r,crfr,c • (33) 

Kotliar and Ruckenstein have shown that, while a cer- 
tain degree of arbitrariness is allowed in the choice of the 
bosonic operator z T , the choice 

Zr,a = (l - d\d r - pl i(T Pr,cr) 2 (elp r .cr + p\_ a d^j 

(l-ele r -pl_ a p r ^y~\ (34) 



turns out to give the correct free fermion limit for.XX — > 
at the mean-field level. In addition, it was shownEa that 
Eq. (|34|) gives the correct non-interacting limit at all 
orders in the fluctuation (1/N) expansion as well. Re- 
striction to the physical Hilbert space is achieved by pro- 
jecting out unphysical states. In a functional-integral 
formalism this is taken care of by integrating over La- 
grange multipliers X^\\r% for each site r@. The total 
effective Hamiltonian is thus given by Eq. ( p7[ ) , whereby 
the Fermi fields have undergone the transformation Eq. 
(|33|), plus the Lagrange multiplier terms: 

H e ff = H t ot + Y i ^4e r + Y.Pr.trPr,a + d\d r - 1^ 

+ Y, 1 X i 2 l (fL.fr., - ptcrPr,* - d\d r ) (35) 
r.cr 



A. Antiferromagnetic phase 

At the lowest order, the functional integral can 
be treated by means of a saddle-point approxima- 
tion, whereby the bosonic fields are replaced by time- 
independent complex numbers. In the antiferromagnetic 
phase, one has the usual two-sublattices saddle point so- 
lution, i.e. 

d r (r) =d, 

e r( T ) = e > 

= »Ai , 

A$ = iX 2 + ^e^ r aA3 , (36) 

from which it follows that 

ZtAt) = z 1 + e^ r az 2 , (37) 

where d, e,pi,p 2 , Ai, A2, A3, zt, z 2 can be chosen to be real 
numbers. 

The effective slave-boson mean-field Hamiltonian be- 
comes 

#SBMF = NX ie 2 - NX 1 + 2N(p 2 +p\){X 1 - X 2 ) 

-4N\ 3 p lP2 + XiNd 2 + NUd 2 - 2NX 2 d 2 

+ Y AU(A2 - M)/k, CT + A 3 Y (T fl+Q,<rfk. t T 
k.tr k,tr 

+ (*? + 4) Y B ( k )/i < r/k,a + H SO(5),I ■ ( 38 ) 

The slave-boson approach can only be applied to deal 
with the on-site repulsion part U , while the SO (5) part 
of the interaction [-^30(5) ji cf. Eq. (|l~8[)1 will still 
be treated by the Hartree-Fock decoupling, as in Eq. 
(|aJ). This is appropriate, since the SO(5) interaction 
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is assumed to be relatively weak (of the order of J) in 
contrast to the large Hubbard repulsion U. In our pro- 
cedure, we first carry out the Hartree-Fock decoupling of 
the SO (5) interaction part of the Hamiltonian, and then 
apply the slave-boson mapping Eq. ( |33| ) . In this way, in 
the AF phase, we have 



H, 



SO(5) 



I — 2-1 k CTC k,<T C k+Q,<7 ' 



with 



,(4) 



v v - 

irjr\ d k\ 2^ Mk'K4,^c k , +Q (T ) 



2 N 



(39) 



(40) 



The transformation, Eq. (p3[), with the saddle-point val- 
ues Eq. ( |36| ) can then be applied to Eq. ( |39| ) by first 
transforming into real space. Going back to momentum 
space yields 

H S0(5)J = E A k ^ {(*? + 4)fLh +Q ,a 
k,<r 

+2z 1 z 2 aflj Ka } (41) 

We can now insert the SO(5) interaction in Eq. j4l| ) 
into the Slave-boson Hamiltonian Eq. fl38|), and obtain 

-H"SBMF = ^k.tr/k.CT 



[e(k)(z 2 - z 2 ) - n + A 2 + 2A| c 4) z 1 z 2 



E /k,(r/k+Q,cr 
k,cr 

AF / 

E r - /k ' 

k,cr 

+Hb os , 



A 3 <t + A k 4 V(z 2 



Hi 



6 OS 



V /k+Q,CT 



?7i(k) ^ (k) 
?/ 2 (k) 773 (k) 



/k,<T 
/k+Q,. 



(42) 



with 



m(k) 

%(k) 
%(k) 



e(k)(z 2 - z 2 ) - n + A 2 + 2A k 4) ziz 2 

A 3 a + A k 4) a(z 2 + z 2 ) 

-e(k)(z 2 - zf) - /x + A 2 + 2A k 4) ziz 2 



(43) 



where Hb os is a purely bosonic part. Diagonalization of 
Eq. ( ft2] ) yields the eigenvalues 

£ ± (fc) = A 2 - /i + 2A< c 4) z 1 z 2 ± {A 2 + e(k) 2 (z 2 - z 2 ) 2 
+2A( 4 )A3(z 2 + z 2 2 ) + (Af) 2 (z 1 2 + z 2 ) 2 } i 

(44) 

Although Eq. (|4|) seems to break the particle-hole 
symmetry at half filling, this is not the case. Indeed 
the self-consistent solution yields at half filling A 2 — fi = 



z 2 = 0, which makes the expression particle-hole sym- 
metric. Eq. (Q) describes the gap structure: Since 
along the magnetic zone boundary the tight-binding dis- 
persion e(k) is identical to zero, the dispersion along the 
magnetic zone boundary reads at half filling: 



£ 2/1 (k) = ±X 3 + A^zf 



(45) 



Therefore, we obtain a constant part of the gap Ajj = 2A 3 
plus a modulation 



A^zf = A mod |d k | , 



with 



^mod 



V 



2 Nz; 



(46) 



(47) 



k.CT 



We introduce the unitary matrix U(k) that diagonal- 
izes the Hamiltonian Eq. (E3): 



(48) 



The self-consistent gap equation is given by Eqs. ( 46 147| ) 
and takes the form (at half filling) 



A 



(4) 



V 



l^klE l d k'| CT (c k ',<r C k'+Q, 



k',<r 
AF 



2 N 



k',<7 



, 2ziz 2 ct z 2 + z| 
U[k) I z 2 + z| 2z!Z 2 a 



The ground-state energy 



*7(k)t 



(49) 



i.i 



AF 



E^ k ) 



(50) 



then has to be minimized with respect to the saddle-point 
values of the fields d, e,px,P2, Ai, A 2 and A 3 . Along with 
the gap-equation Eq. (|4|) this gives 8 coupled equations 
that are readily solved numerically. 



B. SC phase 

The Slave-boson calculation in the superconducting 
phase is very similar to (in fact simpler than) the one 
in the AF phase. Therefore, we will skip details and only 
show the main differences. Indeed, in this case we can 
use the paramagnetic ansatz for the bosons which avoids 
the complication of two different sublattices: 
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d r (r) 
e r (r) 

VrA T ) 
A« 
A (2) 



d 
e 
V 

i\\ 
i\ 2 



gap SC phase (slave boson results) <n>=0.84 



(51) 



z r is obviously independent of a and is given again by Eq. 
(|34|). The saddle-point approximation yields the mean- 
field Hamiltonian (the counterpart of equation (|38"1)): 



#sbmf = A^Aie 2 - TVAi + 2iVp 2 (Ax - A 2 ) 
+ 7Vd 2 (Ai - 2A 2 + f7) 

+ ]T/ 1 L(e(k)z 2 + A 2 - M )/ kiCT 




SO(5),/' 



FIG. 7. The d-wave gap Asc in the superconducting phase 
with density (n) = 0.84 as a function of the Hubbard interac- 
''tion U for fixed V = 0.61 t. We see that, in contrast to the 
antiferromagnetic phase, U lowers the value of the gap. 



In the SC phase, the SO(5) part reads 



H. 



SO(5), 



k.l 



(53) 



with the condition 
A k } = -^7 d k^rfk'(4',T c -kM - <V,T c -k',l) • (54) 



The purely fermionic part 

^F = ^/L[e(k)z 2 + A 2 - M ]/ k ^ 

k,(r 

' ;2 X A k (/ k , T /ik > ;-/k ) T/-k,i 



(55) 



can be diagonalized by the usual Bogoliubov transforma- 
tion. This yields the eigenvalues 



± E(k) = ±\/(e(k)z 2 + A 2 - M ) 2 + (A k 1} ) 2 . (56) 



Calculating the expectation value of the order parame- 
ter in the ground state, constructed from the Bogoliubov 
eigenstates, gives a self-consistent gap equation in the 
usual form: 



i = — z 4 Vd; 

9 /V ' 



2 iV 



1 



V(£(k)z 2 + A 2 - ^ + A^? 



where Asc is the amplitude of the modulation 



A k 1} = d k A sc ■ 



(57) 



(58) 



The gaps obtained from the slave-boson calculation are 
reported in Figs. and ^[ 



Figure g shows the behavior of the AF gap modulation 
as a function of U. Just like in the simp le Hartree-Fock 
type of mean field described in Sec. A_2, we observe that 
increasing the Hubbard-interaction U enhances the d— 
wave-like modulation A moc j of the AF gap at half filling, 
and, eventually, saturates around U/t » 8 — 10. The dis- 
persion of the gap is of course identical to the Hartree- 
Fock mean-field result, Fig. [|. 

The superconducting phase has again a <i-wave struc- 
ture. However, in contrast to the simple Hartree-Fock 
mean field, where U was of no influence, the Hubbard in- 
teraction U results in a suppression of the superconduct- 
ing gap Asc at finite doping, as shown in Fig. ^. This is 
due to the well-known reduction of the effective hopping 
produced by the slave-boson formulation. The combina- 
tion of these two effects with increasing U (enhancement 
of A mo( j in the AF phase and suppression of Asc m the 
SC phase) yields a ratio of the two gaps A mo( j/Asc of the 
order of 10, slightly larger than the Hartree-Fock calcula- 
tion, and more in agreement with experimental findings. 



V. SUMMARY 

Summarizing, we have shown that the |d|-wave-like 
modulation observed in the AF gap at half filling in 
Ca 2 Cu0 2 Cl 2 , similar to the dispersion of the SC gap at 
finite doping, indicates an intimate relationship between 
the two phases, as suggested by the SO(5) theory of high- 
T c superconductivity. The idea is that one gap can be 
mapped into the other by a SO (5) transformation. This 
is done by using an effective SO (5) -invariant Hamilto- 
nian which allows for such a mapping. However, we have 
also shown that it is important to break this symmetry 
by introducing "by hand" an Hubbard interaction term 
in order to correctly obtain the constant part of the AF 
gap. Via this "projected" SO(5) theory, we can interpret 
the experimentally observed d-wave-like modulation of 
the AF gap as the fingerprint of the superconductor in 
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the AF state just like the neutron resonance modeE2r 
may be viewecL-as the fingerprint of the AF correlation 
in the SC statec3. 

The projection can also explain the experimentally 
observed order of magnitude difference between the d— 
wave-like gap modulation in the AF state (w J ~ 
0.12eV) and the d-wave gap in the SC state (ss J/10). 
This different magnitude of the two gap modulations also 
suggests a reason why the Neel temperature differs so 
much from the superconducting transition temperature. 

Finally, we have shown that the projected SO (5) the- 
ory pan provide a qualitative explanation for the observa- 
tion!^, that the maximum T c observed in a large variety 
of high-T c cuprates scales with the next- nearest- neighbor 
hopping matrix element t' . 
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APPENDIX A: DETAILS OF THE MEAN-FIELD 
CALCULATIONS 

1. BCS solution of the SC state 

At finite doping, the symmetry between SC and AF is 
broken to favor SC. Therefore, we look for a BCS solu- 
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tion given by a finite mean-field expectation value of the 
operator A in Eq. (|l^). The Hamiltonian, Eq. (17), 
becomes quadratic in the fermionic operators and can 
be solved by the usual Bogoliubov transformation. The 
mean-field condition results in the BCS gap equation for 
the amplitude of the order parameter Age, i- e. 



H MF = E £ ( k ) c L c k,<x - ^E< nk '>] E 



A 2 d 2 
^SC"k 



(Al) 



which has to be solved in a self-consistent way for a given 
value of the parameter V. Notice that the on-site inter- 
action U does not enter this equation, due to the fact 
that, in contrast to s-wave, two holes of a d-wave pair 
never occupy the same site. 

The fc-dependent SC gap is given by Ascdk and de- 
pends on two parameters: the magnitude Asc and the 
parameter & introduced in Eq. ([^, giving its shape. Our 
procedure consists in fitting these two parameters to the 
recent data on the SC gap by Mesot et al.o. In that 
paper, the SC gap obtained by ARPES experiments was 
fitted with a d-wave form including nearest- and next- 
nearest-neighbor terms: 

rfk = b(cosk x — cos k y ) + (1 — 6)(cos2fc x — cos2fc a ) . 

(A2) 

This form cannot be made consistent with our SO (5) hy- 
pothesis, since it would require a form factor Eq. rtq) 
without the property <?k+Q = — 5-k- However, a fit of 
the SC gap by the form, Eq. (||), turns out to be as 
good as one from the next-nearest-neighbor form, Eq. 

Our best fit to the data of Ref. |20| with Eq. (|9j) gives 
b = 0.81, and Asc = 0-04 f. T he values of the parameters 
can be inserted in Eq. (|Al[) to obtain the appropriate 
value of V (V = 0.89i), which can now be used to study 
the AF phase. 



2. Hartree-Fock solution for the Antiferromagnetic 
phase 

We now look for the AF solution of the mean-field 
equation for the Hamiltonian, Eq. (|i~7|), at half filling. 
This is s lightly more complicated than the BCS result, 
Eq. (Al), due to the interplay between the SO (5) inter- 
action in Eq. ( |l3| ) and the Hubbard repulsion term. 

We introduce a momentum dependent SDW operator 
polarized in the z-direction: 



n k = = 2AT 3 (k) 



(A3) 



Within a mean-field decoupling, the Hamiltonian in Eq. 
(Il7|) becomes 



k.cr 



2 N 



»k 



^LEK'IK<)]5>k 



2 N 



k' k 

E^L/w-a^ 



»k 



«k 



-A mod \d k \n k , 



with the self-consistent parameters 
U 



—[Y< 



2N 
V 

2N 



k' (nv 



(A4) 

(A5) 
(A6) 



It is convenient to recast equation (A4) in matrix 
form as a sum over the antiferromagnetic Brillouin zone 
(AFBZ), 



AFBZ , N f 

Ck.cr ' 



k,<r 



C k+Q,cr 



M ( k).( ck c -J, (A7) 



with 
M(k) = 



e(k) -cr(A t/ +A mod |d k | 
-a(Au + A mod |d k |) -e(k) 



(A8) 



The Hamiltonian (equation ( A7)) can be diagonalized by 
the usual transformation into Bogoliubov operators 



7k, CT = Uk c k.o- + CVkCk+Q,<r 
7k,(7 = '"kCk.cr - <ru k c k+ Q tCr 



(A9) 



with amplitudes 



»k = 



"k = 




(A10) 

and eigenvalues 

± E(k) = ± Ve(k) 2 + (Ay + |d k | A mod )2 . (All) 

At half-filling, the ground state consists of all states of 
7 C particles occupied (c stands for conduction band with 
E(k) < and v for valence band with E(k) > 0) and the 
expectation value (n k ) can be readily evaluated. 
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